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We discuss several examples of how the transverse spatial distribution of partons in the nucleon, as well 
as multiparton correlations, can be probed by observing hard processes (dijets) in high-energy pp (pp) 
and pA (dA) collisions. Such studies can complement the information gained from measurements of 
hard exclusive processes in ep scattering. The transverse spatial distribution of partons determines the 
distribution over pp impact parameters of events with hard dijet production. Correlations in the transverse 
positions of partons can be studied in multiple dijet production. We find that the correlation cross section 
measured by the CDF Collaboration, (T e ff = 14.5 ± l.Tta^ can rje explained by "constituent quark" 
type quark-gluon correlations with r q ~ rjv/3, as suggested by the instanton liquid model of the QCD 
vacuum. Longitudinal and transverse multiparton correlations can be separated in a model-independent 
way by comparing multiple dijet production in pp and pA collisions. Finally, we estimate the cross section 
for exclusive diffractive Higgs production in pp collisions at LHC (rapidity gap survival probability), by 
combining the impact parameter distribution implied by the hard partonic process with information about 
soft interactions gained in pp elastic scattering. 



1 Introduction 

Hard exclusive processes in ep scattering allow to probe not only the distribution of partons with respect 
to longitudinal momentum, but also their spatial distribution in the transverse plane. Examples include the 
hard electroproduction of light mesons and real photons (deeply virtual Compton scattering), as well as the 
photoproduction of heavy quarkonia (J/ip, ip', T). Thanks to QCD factorization theorems the amplitudes 
for these processes can be separated into a "hard" part, calculable in perturbative QCD, and "soft" parts 
characterizing the non-perturbative structure of the involved hadrons. The information about the nucleon is 
contained in so-called generalized parton distributions (GPD's). These are functions of the parton momen- 
tum fractions, x and x', as well as of the invariant momentum transfer to the nucleon, t, and thus combine 
aspects of the usual parton distributions, measured in inclusive deep-inelastic scattering, with those of the 
elastic nucleon form factors. For x = x', their Fourier transform with respect to t describes the spatial 
distribution of partons in the transverse plane. The GPD's thus, in a sense, provide us with a "3D parton 
image" of the nucleon. 

Exclusive processes in ep scattering, however, are not the only reactions which probe the "3D" parton 
distributions. In fact, a lot more information about the longitudinal momentum and transverse spatial 
distribution of partons, as well as about multiparton correlations, can be obtained from the study of selected 
hard processes in (not necessarily exclusive) pp and pA scattering. Comparative studies of ep and pp/pA 



adp data will be provided by the publisher 



2 



Frankfurt, Strikman, and Weiss: 3D partem imaging 



induced hard processes will help to improve the quantitative description of both classes of processes and 
offer many new, fascinating insights into the partonic structure of the nucleon. The conceptual basis for 
such a program is the combination of Gribov's space-time picture of hadron interactions at high energies, 
which allows for a unified description of ep and pp scattering 1 1 1, and the QCD factorization theorems for 
hard exclusive processes. In view of the planned new experiments in both ep scattering (Jefferson Lab at 
12 GeV, EIC/eRHIC) and pp/pA scattering (RHIC, LHC), such studies are very timely. 

In this paper we study various examples of hard processes in pp and pA scattering which probe the 
"3D distribution" of partons and their correlations in the nucleon. First, we show how the production of 
multiple dijets in pp events can resolve spatial correlations of partons in the transverse plane. These cor- 
relations provide interesting new information about the long-wavelength structure of the nucleon, which 
can e.g. be related to the role of instanton-type vacuum fluctuations in generating constituent quark masses 
(spontaneous breaking of chiral symmetry), which is the basis for the chiral quark-soliton model of the 
nucleon 1 8 1. Second, we discuss how similar experiments with nuclear targets (pA, dA) can help to sepa- 
rate between longitudinal and transverse parton-parton correlations in the nucleon wave function. Third, 
we discuss the role of the transverse spatial distribution of partons in the diffractive production of heavy 
particles (Higgs bosons) in pp collisions at LHC energies. Such diffractive events involve a delicate inter- 
play of "hard" and "soft" processes (survival of the rapidity gaps), and the latter can be modeled relying 
on information gained from pp elastic scattering. 



2 The transverse spatial distribution of partons and hard processes in pp 
scattering 

In order to define the transverse spatial distribution of gluons in the nucleon, it is convenient to write the 
gluon GPD in the form 1 (analogous expressions apply to the quark flavor singlet and non-singlet distribu- 
tions) 

H g (x,t,Q 2 ) = g(x,Q 2 ) F g (x,t,Q 2 ), (1) 

where g(x,Q 2 ) = H g (x,t — 0,Q 2 ) is the usual gluon density, and F g (x,t,Q 2 ) the "two-gluon form 
factor" of the nucleon, F g (x,t = Q,Q 2 ) = 1. In the following, the Q 2 -dependence of the form factor 
and related quantities will not be indicated explicitly. One can represent this form factor as the Fourier 
transform of a function of a transverse coordinate variable, p, 

F g (x,t) = J d 2 pe i(Aip) F g (x,p), (* = -Ai). (2) 

cf. the well-known relation between the nucleon electric form factor and the charge density in the Breit 
frame. The function F g (x,p) describes the spatial distribution of gluons with longitudinal momentum 
fraction x in the transverse plane, with J d 2 pF g (x, p) = 1. A measure of the gluonic transverse size of 
the nucleon for given x is the average 



7<~: = / 'rrn 1 F„(.r.p) = 4^- Fg (x,t) t=0 . (3) 



On general grounds, the gluonic transverse size is expected to grow with decreasing x. Different phys- 
ical mechanisms are responsible for this growth in different regions of x. Near x = 1, the growth of (p 2 ) 
is governed by the Feynman mechanism — the t-dependence of the two-gluon form factor disappears if 
the active gluon carries the entire longitudinal momentum of the nucleon, x—*\, causing (p 2 ) to vanishes 
at a; = 1. Note that this behavior is a trivial consequence of the relativistic kinematics of the form factor 



1 We are considering here the GPD for zero longitudinal momentum transfer; in the general case the GPD would depend on the 
momentum fractions of the two gluons separately, x' ^ x. 
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calculation; it does not reflect the transverse position of the x — ► 1 parton relative to the transverse re- 
gion occupied by the spectator system. Hence, in a proton-proton collision (see below) there is no simple 
relation between F g (x, p) and the pp impact parameter in the x-region dominated by the Feynman mecha- 
nism. When x is decreased below the valence region, a distinctive increase of (p 2 ) is caused by pion cloud 
contributions to the gluon density, which set in for x < M^/Mm j2||3]. Finally, when x is decreased 
further, the transverse size grows due to the random walk character of successive emissions in the partonic 
ladder (Gribov diffusion) IQ. 2 

As to the Q 2 dependence, the nucleon's gluonic transverse size at fixed x decreases with increasing Q 2 
as a result of DGLAP evolution; see Ref. |4| for a discussion of this effect. 

In ep scattering, the cleanest way to probe the transverse spatial distribution of gluons is via exclusive 
photo- or electroproduction of heavy quarkonia (J/tp, T). These processes have been measured in a num- 
ber of fixed-target experiments (x > 10 _1 ), as well as in the HI and ZEUS experiments at the HERA 
collider (x ~ 10~ 2 — 10~ 3 ), see Ref. |5| for a recent review of the data. The two-gluon form factor can be 
extracted from the measured t-dependence of the differential cross section. It was found that (p 2 ) increases 
from ~ 0.24 fm 2 at x ~ 10 _1 to ~ 0.35 fm 2 at x ~ 10~ 2 — 10~ 3 ; this difference can quantitatively be 
explained by the pion cloud contribution which is switched on for x < M^/Mm Q. Concerning the shape 
of the two-gluon form factor, it has been argued that at x > 10 _1 (where pion cloud contributions are ab- 
sent) the two-gluon form factor should follow the axial form factor of the nucleon, and thus be described 
by the dipole parametrization 

F g {x,t) = {l-t/m 2 g )- 2 , m 2 = 1.1 GeV 2 (x > 10" 1 ). (4) 

This form indeed describes well the i-dependence of the fixed-target J/ip photoproduction experiments. 
Based on Eq. 10), we have suggested in Ref. |4| a generalized dipole parametrization valid also at small x, 
which incorporates the observed increase in the gluonic transverse size (as well as the effects of DGLAP 
evolution) by way of an x— and Q 2 -dependent of the dipole mass parameter, m 2 (x,Q 2 ). This simple 
parametrization conveniently summarizes our knowledge of the transverse spatial distribution of gluons 
from ep scattering. For details, see Refs. |4]|5]. 

The transverse spatial distribution of quarks in the nucleon can be probed in the production of neutral 
vector mesons at sufficiently large x, and also in processes where gluon exchange does not contribute, such 
as the production of charged vector mesons like 7* + p — > p + + n. 

Turning now to pp collisions, an immediate application of the transverse spatial distribution of partons 
is in the description of the impact parameter dependence of the cross section for hard dijet production 
|4|. In a pp collision with cm. energy y/s, a hard dijet with transverse momentum q± at zero rapidity is 
produced in the collision of two partons carrying momentum fractions 

Xi =x 2 = 2q ± /y/s (5) 

of the respective protons; the generalization to non-zero rapidity, Xi 7^ x 2 , is straightforward. The proba- 
bility for such a parton-parton collision, as function of the impact parameter of the pp system, b, is given 
by the convolution of the spatial distributions of the partons; for gluons 

Pi{b) = J d 2 Pl J d 2 P2 5^(b- pi+p 2 ) F g (x uPl ) F g {x 2 ,p 2 ). (6) 

The scale of the parton distributions here is qj_. One can easily generalize this to the production of multiple 
dijets. Neglecting possible correlations in the transverse positions of two partons in the same proton (this 
question will be discussed in detail in Section^, the ^-dependent probability for the production of a double 



For a cross sections increasing as a power of the energy, unitarity at small impact parameters also leads to an increase of the 
radius of the interaction with energy. 
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dijet would be given by 

P 4 (6) = - [P2(6)]2 2 . (7) 

While exclusive processes in ep scattering provide in principle the cleanest way to access the transverse 
spatial distribution of partons, there are several instances in which pp scattering is more effective. One 
is the study of large x, where the cross sections for exclusive processes in ep are small. Besides, in this 
kinematics, when x ~ 1 are probed, one is mostly sensitive to the Feynman mechanism; see the discussion 
above. In this case the global transverse distribution of matter can be measured more directly using var- 
ious reactions combining a soft and hard trigger, in particular in connection with pA collisions |6|. New 
opportunities for such studies will emerge at LHC, where the high luminosity will allow, for example, to 
compare the characteristics of W + and W~ production at the same forward rapidities, corresponding to 
relatively high x where the d/u ratio deviates strongly from the naive value 1/2. By studying the accom- 
panying production of hadrons one can learn which configurations in the nucleon have larger transverse 
size — those with a leading u-quark or with a leading d-quark. One suitable observable is, for example, 
the distribution of the number of events over the number of the produced soft particles. A larger transverse 
size corresponds to a larger probability of soft interactions, and hence to a larger probability of events with 
large multiplicity. It is interesting to note that the studies of the associated soft hadron multiplicity in the 
production of W and Z bosons in pp collisions by the CDF collaboration at Fermilab find an increase 
of this multiplicity by a factor of two as compared to generic inelastic events |7|. This appears natural if 
one takes into account that the hard quarks producing the weak bosons have a narrower transverse spatial 
distribution than the soft partons. As a result, the average impact parameters in events with weak boson 
production are much smaller than in generic inelastic collisions, leading to an enhancement of multiple 
soft and semi-hard interactions |4|. 



3 Probing correlations in the proton parton wave function via multiple di- 
jet production 

Single parton densities and GPDs do not carry information about longitudinal and transverse correlations 
of partons in the hadron wave function. Such information can be extracted from high energy pp and pA 
collisions where two (or more) pairs of partons can collide to produce multiple dijets, with a kinematics 
distinguishable from those produced in 2 — > 4 parton processes. Since the momentum scale of the hard 
interaction, p t , corresponds to much smaller transverse distances in coordinate space than the hadronic 
radius, in a double parton collision the two interaction regions are well separated in transverse space. 
Experimentally, one measures the ratio 

ffl^W ^ + P ~+ jCtl + jCt2 + j6t3 + 7) 



da _.,.„.. da . 

(p + p-»]etl +jet2) • (p + p^jet3 + 7) 



dQidfl2 dQ,zdQ,i 

f(x 1 ,X 3 ,/J, 2 )f(x2,X4,H 2 ) 



(8) 



CToff f(xi,Ll 2 )f(x 2l fl 2 )f(x3, Ai 2 )/(x 4 , fJ. 2 ) ' 

where f(xi,xs), f(x2, X4) are the longitudinal light-cone double parton densities at the hard scale [x 2 (we 
assume for simplicity that the virtuality in both hard processes is comparable; in the following equations 
we suppress dependence on /i 2 ), and the quantity a c g can be interpreted as the "transverse correlation 
area". The variables £1; characterize the observed jets (or photons) — their transverse momenta, rapidities, 
cuts on the opening angle, etc. 

Parton correlations can emerge due to nonperturbative effects at a low resolution scale, or due to the 
effects of QCD evolution. One possible nonperturbative mechanism is the existence of "constituent quarks" 
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within the nucleon, which appear due to the interaction of current quarks with localized non-perturbative 
gluon fields, resulting in local short-range correlations in the transverse spatial distribution of gluons. The 
instanton model of the QCD vacuum suggests a constituent quark radius of about 1/3 the nucleon radius, 
r q ps tn/3. Another nonperturbative mechanism, relevant at small x, are fluctuations of the color field in 
the nucleon due to the fluctuations of the transverse size of the quark distribution. Perturbative correlations 
emerge due to small transverse distances in the emission process in the perturbative partonic ladder in 
DGLAP evolution. Of all the mentioned mechanisms, only the first one is effective at x > 0.05, where the 
data of the CDF experiment were collected. 

The CDF experiment observed correlation effects in a restricted x range (two balanced jets, and jet plus 
photon) and found a e s = 14.5 ± 1.7^2 3 m ^>- This value is significantly smaller than the naive estimate 
obtained by taking a uniform distribution of partons of a transverse size determined by the e.m. form factor 
of the nucleon, which gives a e g as 53 mb, indicating strong correlations between the transverse positions 
of partons in the transverse plane. The longitudinal correlation between partons in the measured kinematics 
due to energy conservation is likely to be small, as x% + X2 and X3 + X4 are much smaller than 1. If this 
effect were important it would likely lead to a suppression of the double parton collision cross section, and 
hence to an increase of <r e ff . However, no dependence of <r e g on Xi was observed in the experiment. 

For a more quantitative analysis of the CDF data, we can make use of the information about the trans- 
verse spatial distribution of gluons gained from J/ip photoproduction, as summarized in Section|2] Since 
the x values of the partons probed were reasonably small compared to 1, the simple "geometric" picture of 
the pp collision in transverse position in the spirit of Eqs. and is justified, and one has 



Ccff 



(9) 



Evaluating this with the dipole parametrization of the two-gluon form factor this comes to 



28tt 

a eS = -— « 34 mb. (10) 

g 

Thus, about 50% of the enhancement compared to the naive estimate of the previous paragraph is due to 
smaller actual transverse radius of the gluon distribution. Still, our value indicates significant correlations 
in the transverse positions of the partons. In the kinematics discussed here the relevant partons are both 
quarks and gluons. We can estimate the effect of correlations assuming that most of the partons are con- 
centrated in a small transverse area associated with the "constituent quarks", as implied by the instanton 
liquid model of Diakonov and Petrov [8 1. Assuming a constituent quark radius of r q ~ rjy/3, we obtain 
an enhancement factor due to transverse spatial correlations of partons of 



1 



- + -^ ~ 1.6-2. (11) 

9 ri 



This is roughly the value needed to explain the remaining discrepancy with the CDF data. Thus, the 
combination of the relatively small transverse size of the distribution of large-x gluons and the quark- 
gluon correlations implied by "constituent quarks" with r q w rjv/3 is sufficient to explain the trend of the 
CDF data. Further studies of multijet events at hadron-hadron colliders, with a broader range of final states, 
would in principle allow to measure separately quark-quark, quark-gluon, and gluon-gluon correlations 
for different x. 

However, studies based on pp or pp collisions alone do not allow for a model-independent separation 
of transverse and longitudinal correlations. This is possible only in pA collisions at RHIC and LHC. The 
reason is that the nucleus, having a thickness which practically does not change on the nucleon transverse 
scale, provides an important contribution which is sensitive only to the longitudinal correlations of hadrons 
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|9|. This is the contribution when two partons of the incident nucleon interact with partons belonging to 
two different nucleons in the nucleus, 02, 



(J-2 



NN 
double" 



A-l 



A 



J d 2 bT 2 {b) 



(12) 



The other term is the impulse approximation — two partons of the incoming nucleon interact with two 
partons of the same nucleon in the nucleus, <7i, which is simply equal to A times the cross section of 
double scattering in pp collisions. Thus, by measuring the ratio of pA and pp double scattering cross 
sections we can determine 1 / er e g via the relation 



1 



pA 
^double 

^"double 



1 - 1/A 



Jd 2 bT 2 (b) 



(13) 



This expression applies for xa > 0.03, where nuclear effects in the structure functions are small. Note 
that the experimental measurement of the A-dependence will provide an independent test of this equation. 
For the ratio of double to single scattering terms we find, for A > 12, 



R = 



1l 
01 



0.68 



0.:SS 



geff 
14 mb ' 



(14) 



Taking the CDF value of a c g ~ 14 mb, we obtain R ~ 3 for A ~ 200. Thus, the separation of the two 
terms will be quite straightforward. Even in the case of deuteron-nucleus scattering, which was studied 
at RHIC recently, the contributions from two partons of one nucleon of the deuteron interacting with 
two different nucleons in the nucleus remains significant.lt constitutes about 50% of the cross section for 
A ~ 200. Hence we conclude that both measurements of pA and dA collisions will allow to measure of 
cr e ff if it is > 5 mb, with pA being a better option. Finally, if <r c ff will have been measured in pA collisions, 
it will be possible to extract the longitudinal two-parton distributions in a model independent way. 

To summarize, we have demonstrated that future experiments will be able to measure independently the 
longitudinal and transverse two-parton distributions in the nucleon. With a detector of sufficiently large 
acceptance it would be possible to extend these studies even to the case of three parton correlations. 



4 The transverse spatial distribution of gluons and diffractive Higgs pro- 
duction at LHC 

Exclusive diffractive production of Higgs bosons, 

p + p -> p + (gap) + H + (gap) + p, (15) 

seems to be one of the promising candidates for the Higgs search at LHC; see Ref. 1 10 1 and references 
therein. From the point of view of strong interactions, such processes involve a delicate interplay between 
"hard" and "soft" interactions, which leads to a characteristic dependence of the cross section on the impact 
parameter of the pp system, b. The heavy particle is produced in a hard partonic process (virtualities 
~ Mjj) involving the exchange of two gluons between the nucleons — one for the gluon-gluon fusion 
making the Higgs, the other for color neutralization. The impact parameter distribution for this process is 
described by the function P4Q)), see Eq. (0. In addition, the soft interactions between the two nucleons 
(viz. the spectator systems) have to conspire in such a way as not to fill the rapidity gaps left open by the 
hard process. The probability for this can be inferred from the amplitude of pp elastic scattering in the 
impact parameter representation, T(s, b). Namely, 



|i-r( Sj 6)| 2 



(16) 



adp header will be provided by the publisher 



7 



0.5 - 




b [fm] 



(a) 



»5 



CM 



0.06 



0.04 



0.02 



Islam et al. 
Khoze et al. 




b [fm] 



(b) 



Fig. 1 (a) The probability distribution for no inelastic interaction, Eq. 1 1 6i . for y/s = 14 TeV, as calculated with 
the parametrizations of the pp elastic amplitude of Islam et al. 1111 . and Khoze et al. 1121 . (b) The fo-distribution 
for diffractive exclusive Higgs production, Eq. 1171 . as obtained with the dipole-type two-gluon form factor with 
Shown are the "radial" distributions including a factor 2-Kb. 



is the probability for having no inelastic interaction in a pp collision with impact parameter b. This function 
can be evaluated using available phenomenological parametrizations, which can be extrapolated to the 
LHC energy, y/s = 14 TeV. Fig.[U shows |1 — T(s,6)| 2 for the parametrization of Ref.|ll| and the 
multipomeron model of Ref. 1 12 1. Both parametrizations indicate that at this energy the nucleon is "black" 
[|r(s, 6)| ss 1] for b < 1 fm. Fig.^3 shows the product 

\l-T(s,b)\ 2 P,(b), (17) 

which governs the 6-distribution of the cross section. The distribution is suppressed at small b (because 
the probability for no inelastic interaction is small) as well as at large b (because the overlap of the gluon 
distributions, Pi{b), vanishes), and is thus concentrated at intermediate values of the impact parameter, 
b - 1 fm. 
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Fig. 2 The rapidity gap survival probability, S 2 , Eq. i 1 SI . obtained by integrating the product |1 — T(s, b)\ 2 Piib) 
shown in Fig.Qover impact parameters. Shown is the result as a function of s, for various values of the mass parameter 
in the two-gluon form factor, m 2 . The Tevatron and LHC energies are marked by arrows. 



The integral of Eq. d!7l > defines the so-called rapidity gap survival probability, 

S 2 = J d 2 b |1 -r(s,6)| 2 P A {b). (18) 

Fig-Hshows the variation of this quantity with s between Tevatron and LHC energies, for various values 
of the dipole mass in the two-gluon form factor of the nucleon, Eq. 0. The gap survival probability 
decreases with s because the "black" region in the proton grows with the collision energy. Note that our 
results for S 2 are in agreement with those obtained by Khoze et al. 1 12 1 in a multi-pomeron model, given 
the uncertainty in the basic nucleon size parameter (denoted by b) in that approach, as well as with those 
reported by Maor et al. 1131 . In view of the different theoretical input to these approaches this is very 
encouraging. It is worth noting that, if we consider the production of an object of fixed mass at different 
energies, the values of Xi decrease with the energy, corresponding to a smaller effective value of m 2 in 
Fig. |2| When this is taken into account, the actual drop of the survival probability with energy become 
much more modest. 
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